Structural phase transitions in low-dimensional ion crystals 
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A chain of singly-charged particles, confined by a harmonic potential, exhibits a sudden transition 
to a zigzag configuration when the radial potential reaches a critical value, depending on the particle 
number. This structural change is a phase transition of second order, whose order parameter is the 
crystal displacement from the chain axis. We study analytically the transition using Landau theory 
and find full agreement with numerical predictions by J. Schiffer [Phys. Rev. Lett. 70, 818 (1993)] 
and Piacente et al [Phys. Rev. B 69, 045324 (2004)]. Our theory allows us to determine analytically 
the system's behaviour at the transition point. 
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INTRODUCTION 



Wigner crystals of ions in Paul or Penning traps are a 
remarkable example of selforganized matter at ultralow 
temperatures [l[ . These systems are usually composed of 
singly-charged particles, which are kept together by ex- 
ternal time-dependent radio-frequency or static magneto- 
electric potentials, and which reach crystallization by 
means of laser cooling. Among several important as- 
pects, the transition from disorder to order for few ions 
was studied in 0, S 3 ; long-range order in three dimen- 
sional structures in Penning traps was first demonstrated 
in 0, 0|; and more complex crystalline structures have 
been realized, see for instance 0, H, 0- Most recently, 
these crystalline structures have been attracting increas- 
ing attention for the realization of quantum information 
processors 



11, 121 and simulators hi 14, 15, h 



In this perspective, the clear understanding and char- 
acterization of the structural properties would provide 
the possibility to control at the microscopic level the 
dynamics of complex systems. Moreover, ion crys- 
tals are systems characterized by truly long-range, un- 
screened Coulomb interactions, and hence constitute in- 
teresting physical systems where one can test equilibrium 
and out-of-equilibrium statistical mechanics models for 
systems exhibiting non-extensive thermodynamic func- 
tions 

Structural transitions in ion crystals are induced ei- 
ther by changing the external potentia l flil . [20j or by 
introducing other forms of instabilities |2l| . Structural 
transitions in low dimensional ion crystals were first char- 
acterized experimentally in [l^, [2^. Here, starting from 
a chain configuration, the sudden transition to a planar 
zigzag structure, as shown in Fig. (TJ was observed when 
the radial potential reached a critical value, dependent 
on the ion number. In theoretical investigations it was 




FIG. 1: Structural phase transition in a string of equidistant 
trapped ions from a linear to a planar zigzag configuration. 
For ions in a harmonic trap, close to the transition point the 
zigzag configuration is evident about the center of the trap, 
where the density of ion is larger [2^ ]. 



conjectured that the structural change from a chain to 
a zigzag is a second order phase transition [2^ . Further 
numerical work showed that at this transition point the 
ground state energy is characterized by a discontinuity in 
the second derivative with respect to the particles den- 
sity 

In this article we study the structural phase transition 
of an ion crystal from a linear chain to a zigzag configu- 
ration in a suitably defined thermodynamic limit, by de- 
veloping an analytic theory which allows us to determine 
the behavior of the system at the critical point. From 
symmetry considerations we conjecture the spontaneous 
symmetry breaking 
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we 



Applying Landau theory 
identify the order parameter and the soft mode driving 
the instability, and demonstrate that the system under- 
goes a second order phase transition. Our theory is valid 
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at T = 0, when the system exhibits long-range order. 
It allows us to determine the system's behavior at the 
transition point, and the results we find are in agreement 
with the numerical results reported in 22 1 and in [ist . 

This article is organized as follows. In Sec. we in- 
troduce the model and discuss first the transition for a 
chain of 3 ions from a linear to a zigzag configuration 
of charges. In Sec. we derive the dispersion relations 
and eigenmodes at equilibrium of the linear chain and 
of the zigzag configuration in the thermodynamic limit. 
In Sec. we focus onto the classical phase transition be- 
tween the two configurations, identify the soft mode and 
study analytically the system around the critical point. 
In Sec. we conclude and in the appendices wc report the 
details of calculations presented in Sec. . 



ORDERED STRUCTURES OF IONS IN LOW 
DIMENSIONS 



The model we consider is constituted by N particles of 
mass m and charge Q, which arc confined by an exter- 
nal harmonic potential along one axis. The particles are 
classical, and the Hamiltonian governing their dynamics 
reads 



H 
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(1) 



where Vj — (xj,yj, Zj) and Pj are the positions and conju- 
gate momenta, with j = 1, . . . , iV. The term V accounts 
for the oscillator's potential and the Coulomb repulsion. 
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Here, the potential is characterized by harmonic confine- 
ment at frequency v and i>t h^ the axial and transverse 
direction, respectively, whereby vt > v for the case we 
are going to study. We denote by a = vt/v the trap 
aspect ratio, such that a > 1. 

At sufficiently low temperatures, the ions localize 
themselves at the equilibrium positions rj''^ which solve 
the coupled equations describing the equilibrium of the 
forces, 
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(3) 



When the transverse frequency vt exceeds a critical value 

(c) 

vl , which depends on the axial trap frequency v and on 
the number of ions, the solutions of Eq. ([3]) are aligned 
along the x-axis, forming a string. Tables of the equi- 
librium positions for string up to 10 ions have been re- 
ported in psi . [26j . An analytical form for the linear 



density of ions along the trap axis at equilibrium was 
determined in [S^] for ^ 1 and using the local den- 
sity approximation. The linear fluctuations about the 
classical ground state of an ion chain in a harmonic 
trap have been analytically studied in [2^, |2^. This 
study identified as well the value of the critical trans- 
verse frequency v'f^ k, SNv/ {iy/log N) , using an expan- 
sion at leading order in 1/ log A^, and by considering only 
nearest-neighbour contributions. Within this approxima- 
tion this value is consistent with numerical results [2^ . 
and is in goo d agreement with previous analytical evalu- 
ations in [30|, which calculated the critical value taking 
into account the long-range interaction between the ions 
but assuming that the particles are equidistant. 

When the transverse frequency is varied, so that < 
, the stable configuration is first a zigzag structure, 
then at smaller values it has an abrupt transition to a he- 
licoidal one, and so on thereby acquiring more complex 
structures [3, 0, HO]. Eventually, for a large number 
of ions and for aspect ratios a sufficiently close to unity 
the structure is expected to take the b.c.c. crystalline 
form [ij. In the following, wc study the transition from 
an ion chain to a zigzag structure for the most simple 
model, namely three ions in a linear Paul trap. This 
system allows us to get some insight into the system, 
before considering the structural transition in the ther- 
modynamic limit in Sec. . 



Structural stability of a three-ion chain 

We consider A^ = 3 ions inside a trap with vt > and 
calculate their equilibrium positions as a function of the 
aspect ratio a ~ vtjv. We restrict for simplicity to two 
dimensions, which we here identify with the x — y plane, 
and rewrite the potential ^ in dimensionless variables 
as 
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where x- = Xi/l, y- = yi/l , l'^ 



^'limv^) and V 
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Vjipmv^'). Throughout this section we drop the prime 
superscript. 

The normal modes frequencies for the linear and the 
zigzag structures are displayed in Fig. [5] as a function of 
a, when a is decreased across the value for which the 
linear chain becomes mechanically unstable. The transi- 
tion to the zigzag configuration takes place at the value 
a = a*, such that the smaller transverse frequency of 
the linear chain vanishes. Wc now study in detail the 
classical equilibrium positions for a > 1. Assuming the 
convention x\ < xi < X3, the symmetry of the trap- 



ping potential imposes x^^"* = 0, x\ 
with x > 0. From Eqs. ([3]) we also find 2/2 
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FIG. 2: (color online) Excitation frequencies of the oscilla- 
tions modes of a three-ion string, in units of the axial trap 
frequency v,as,& function of the trap anisotropy a. The eigen- 
modes for a> a* are the longitunal (top) and the transverse 
ones (bottom). For a < a* the eigenmodes are combination 
of longitudinal and transverse modes with opposite parity by 
reflection about s = 0. 



yf'^ = y"^'^ = ^, with y > 0. The hnear configuration, 
namely the set of solutions with y = 0, is found when 
the aspect ratio a > a*, where a* = yjl2/b. When 
a < a*, then y > and the structure becomes planar. 
We denote this case by "zigzag configuration", as it is 
indeed the most elementary instance of the structure one 
observes for many ions. Here, for a < a* terms x, y take 
the form 
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Their functional dependence on the aspect ratio a is dis- 
played in Fig. [3l One can observe the discontinuity of 
the derivative at a*, corresponding to the transition to 
a different equilibrium configuration. For a ^ a* the 
change is faster for the transverse displacement, as it is 
visible by the expansion oi x,y at 5a ^ a* — a, 

y ^ yo da^ + O (^Sa^^^'j 

X = xiin ^ xq 5a + ((5a^) 

where x\\a = (5/4)^/^ is the value taken by x when the 
linear chain is stable, while yg « 0.74, xo « 1.85. We 
note that, about the instability point of the linear chain, 
the transverse displacement y plays the role of the order 
parameter, while the changes of the axial distance x are 
induced by the changes of y, and therefore about the 
value a* these are less dramatic. 




FIG. 3: Equilibrium position of the external ions of a string 
of 3 particles as a function of the trap anisotropy a. The 
solid and dashed lines display the longitudinal and transverse 
variables, x and y, in units of the characteristic length I. The 
vertical dotted line indicates the transition value a*, where 
the equilibrium configuration makes an abrupt change from a 
linear chain to a zigzag structure. 



THE LINEAR AND THE ZIGZAG STRUCTURES 

In this section, we study the static properties of the 
linear chain and of the zigzag configuration in the ther- 
modynamic limit. For an ion chain inside a trap, a good 
thermodynamic limit is found by fixing the interparticle 
spacing a at the chain center when N oo. This corre- 
sponds to the requirement that the axial trap frequ ency 
vanishes according to the relation v ~ v^IoglV /N |28l.l23|. 
In this limit, the critical transverse frequency v^'^^ is con- 
stant, and the behaviour at the mechanical instability is 
equivalent to that of a uniform chain with equal inter- 
particle distance a between neighbouring ions [sij. The 
uniform chain is the model we will use for determining 
the ground state and the motion of the linear and zigzag 
structure in the thermodynamic limit. 



The linear chain 



We assume a stable linear chain of ions, namely vt > 



In this limit the equilibrium positions lie along the 
0, 0), and we use the convention Xi > 



(0) , (0) 

x-axis, r^- — [Xj 
Xj for i > j. For small vibrations around these points we 
approximate the potential in Eq. ^ by its second order 



Taylor expansion in the displacements qj ~ x 
Zj. In this limit the equations of motion are 



2 ^ m 



(6) 
(7) 
(8) 
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FIG. 4: (color online) Excitation spectrum of the uniform 
chain. The eigenfrequencies ui, in units of ujq = jmo? , 
are plotted as a function of the quasimomentum fc, in units 
of 7r/a. The axial spectrum (green solid line) and the trans- 
verse spectrum (red dashed line) are obtained from Eqs. (|lip 



and (|12p. respectively. Here, Vt = 1.1- 



and describe a system of coupled oscillators, with 
long range interaction and position-dependent coupling 



strength, 
read 



Here, the coefficients /C. 
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We note that at second order in the harmonic expansion 
the axial and transverse vibrations are decoupled. It is 
easily verified that the center-of-mass motion is an eigen- 
mode of the secular equations ©-([S]) at eigenfrequencies 
v and Vf for the axial and transverse motion, respectively. 
The solution to Eqs. ©-([S]) have been studied in [2^.i29t. 

For the purpose of studying the behaviour at the me- 
chanical instability, we now consider the simplified model 
of the uniform chain, where the interparticle distance at 
equilibrium is fixed. This case is found by setting i/ = in 
Eqs. ([B])-® and assuming constant interparticle distance 



(0) 



' . Such condition can be realized for the 
central ions of a long ion chain inside of a linear Paul trap 
[3^ or for ions confined in a ring of large radius 3 2o| . 
This second scenario corresponds to take periodic bound- 
ary conditions, qi = Qn+i, etc. Crystallization is found 
assuming, for instance, that one ion is pinned at the po- 
0. Then, the classical equilibrium positions 
are xy = ja and the coupling strengths in Eqs. (Ull-dS]) 



sition X. 
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take the form 



form 
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The dispersion relations are [3S 
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with k = 2-KnlNa and n = 0, ±1, ±2, . . . , iV/2. The 
spectrum corresponding to Eqs. ((TT|) and p2|) is shown 
in Fig. 131 The axial eigenmodes at frequency u\\ (fc) are 



such that 



yii:(ei- 

fc>0 



COS kja -I- 8^ sin kja 
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where the superscript ± indicates parity by reflection 
k —k. Analogously, we denote the transverse eigen- 
modes at frequency w±(fc) by ^^^^^ and where 



Vi 



fc>0 



COS fcja + ^ sin fcja ) ,(14) 



cos fcja -I- 'I'^.^ "^sinfcja) .(15) 



fc>o 



We note that the modes at fc = 7r/2 are even. A close in- 
spection to Eq. (fT2|) shows that there may exist values of 
the transverse trap frequency, at fixed interparticle dis- 
tance a, for which u]\ < 0, that is, imaginary frequency 
solutions. For such values, thus, the chain is unstable. 



The threshold value vf* , such that for vt > v"^' the lin- 
ear chain is stable, is found by solving minfc(wj_) = 
(see Sec. ). The minimum is found at k = 7r/a and cor- 
respondingly 



2Q2 
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(16) 



where result pB]) is found for N ^ 00 using J^eyoC^^ ~ 
= (1 — 2^P) C(p), with ({p) the Riemann-zeta func- 
tion. The value in Eq. (|16|) depends on the interparti- 
cle spacing a and provides the range of validity of the 
results presented in this section. It coincides with the 
value reported in [s^, where a similar model to the one 
discussed here was considered. It is close to the result 
^(c),trap2 _ 4n2/,7^Q^^^p(o)3 found at leading order in 

l/\ogN in [29|, where a{x) gives the interparticle dis- 
tance as a function of x in the local density approxi- 
mation, and a(0) = a is the value at the chain center. 
This result was obtained by considering the inhomoge- 
neous distribution of ions along the chain, but keeping 
only the nearest-neighbours interaction. The small dis- 
crepancy between the two values is to be attributed to 
the different approximations that have been applied in 
each model. 



The zigzag structure 



For h't < vf^ ^ and sufficiently close to the critical value, 
the stable configuration is a lAg'L&'g structure. We now 
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FIG. 5: Transverse equilibrium displacement 6, in units of 
the interparticle spacing a, as a function of the transverse 

(c) 

frequency vt in units of vl . On the right of the curve the 
ion crystal is a linear chain. In the region on the left of the 
curve it exhibits a zigzag structure. 



with X = ^/fl- The cocfScients Eq. ([10]), and the cor- 
responding equations of motion for the linear chain, 
Eqs. (0,(17]), are recovered for x ^ 0. 

In general, the structural change brings to a doubling 
of the unit cell d of the crystal, which from d = a in 
the linear chain goes to d = 2a in the zigzag configura- 
tion. Correspondingly, the Brillouin zone of the zigzag 
is reduced by a factor 2, and the wave vectors now take 
the values k = 2im/Na and n = 0, ±1, ±2, . . . , In 
Eqs. ([T5)) and ([TO]) one can easily verify that the bulk 
excitations are eigenmodes of the chain at frequencies v 
and Vf The other eigenvalues and eigcnfunctions can be 
found using the ansatz fn''^\ with 



evaluate its dispersion relation and eigenmodes for ions 
on a ring and for periodic boundary conditions. We 
assume the equilibrium positions to lie on the x ~ y 



plane with 



(0) _ ^^(0) „(0) 



(xr,2;r,0). Then, x. 



(0) _ 



and 



— (— 1)"6/2, with b a real and positive constant. 



which is determined from the equation 
■nw} v-^ 4 



^[(2£-l)2a2+62]3/2 
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Figure [5] displays the transverse equilibrium displacement 
6 as a function of the transverse frequency vt, as it is 
obtained by solving numerically Eq. (|17[) . Assuming that 
the zigzag configuration is stable, wc denote by qn ~ 
Xn — , Wn = Un — Vn^ i and z„ the axial and transverse 
displacements, and expand the potential of Eq. ^ up 
to second order. In this limit the motion along the z 
direction is decoupled from the vibrations on the plane, 
and the resulting equations of motion for qn and Wn read 

mqn = -^IC'fiqn - q7i+e) (18) 
-{-^r^yeiwn - Wn+e) , 

mWn = -mV^Wn + ^ IC^{Wn - Wn+l) (19) 

-{-irY.ye{qn-qn+,) . 

The coefficients appearing in these equations depend only 
on the interparticle distance, ia = (n' — n)a, as the struc- 
ture is periodic along x. In particular, for £ even they 
read 

2Q2 1 
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while for (. odd they are given by 
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where j ~ 1,2 and ka varies on the interval [— 7r/2, 7r/2]. 
In particular, we note the relation in' ■* (k) ~ fr/'^'' (fc -|- 
ir/a). The corresponding eigenmodes are given by the 
real and imaginary parts of these vectors. Using this 
ansatz, we obtain the coupled equations 



u;,Mkr=c[^\k) + e^^-'^^B{k) 
{v^,^u,,^{kr)^d^\k) 



(21) 

^B{k) , (22) 



whereby 



B{k) = — V 3^2^-1 sin(2£ - l)A:c 
ra ^-^ 
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C^+> [k) = —J2( sin^ kia + /C^^_^ cos^ 
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The eigenfrequencies are found by eliminating the param- 
eter e^ '^^ from Eqs. ((2T|) - (f22]) . The excitation spectrum 
exhibits four branches in the new Brillouin zone, and 
their functional dependence on k is 
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FIG. 6: (color online) Branches of the excitation spectrum 
of a zigzag structure for the modes on the x — y plane, as 
obtained from Eq. (|23|l . The curves display the frequencies 
uj2,+ {k) (green solid), uj2,-{k) (blue dotted), uji^+{k) (orange 
dot-dashed), uji^-{k) (red dashed), in units of ujo, as a func- 
tion of k, in units of tt/o. The Brillouin zone is now half the 
Brillouin zone of the linear chain due to the doubling of the 
crystal periodicity. Here, ut = O.QVf 



with j = 1,2. The spectrum for the excitations on the 
X — y plane is displayed in Fig. [51 

We note that in the limit 6^0 the branches of the 
spectrum of the linear chain, Eqs. and p2|) . are re- 
covered from Eqs. ((23l) . In fact, for 6 = we have B = 



and C[ 



(±) 



2C. 



(T) 



such that each solution has double 



degeneracy, with 



b=0 



b=0 



and 



(fc)^ = (fc)^ = j/f - (fc) /2 



6=0 



6=0 



which reproduce respectively Eqs. ([TT|) and (|T2|) (note 
that 1/2 - '(A:) - C^^^(fc) < 0). 



LANDAU THEORY OF THE STRUCTURAL 
PHASE TRANSITION 

If the ions are crystallized along a line, by lowering 
the transverse confinement vt the system will be led to 
a situation in which the linear chain gets unstable. In 
this regime, one observes experimentally a transition, in 
which the ions are crystallized on a plane, according to 
a zigzag distribution of particles. In the literature it was 
conjectured that this is a second-order phase transition. 
This comecture is supported by the numerical results 



Indeed, one can observe that the transition from a lin- 
ear to a zigzag configuration is characterized by a sym- 
metry breaking resulting in the increase of the unit cell 
by a factor of 2. It is combined with a transition from 
a linear to a planar structure corresponding to the loss 
of rotational symmetry about the x-axis. Then, one can 
identify the order parameter with the displacement of 
the equilibrium position from the x-axis, while the con- 
trol parameter can be taken as the transverse frequency 
i>t when the interparticle distance is fixed. Starting from 
this educated guess we apply Landau theory to the tran- 
sition |2J]. We focus on the situation in which the inter- 
particle distance a is fixed, and study the crystal struc- 
ture when the transverse confinement i>t varies across the 

(c) 

critical value . We explicitly determine the critical ex- 
ponent of the order parameter around the critical value, 
and find that it is in agreement with the numerical results 
in [H. 



The soft mode 

Let us now go back to the dispersion relation for the 
transverse modes of the linear chain in Eq. ((T2|). The 

(c) 

structural transition takes place for the critical value i'^ , 
Eq. ([TO)) , such that the frequency of the lowest transverse 
mode of the linear chain vanishes, as shown in Fig. [7^. 
The smallest transverse frequency lu± is found at the 
value of the wave vector k, at which the semipositive- 
definite function 



N 



r 



is maximum in the interval < if < 7r/2, as seen from 
Eq. ini) for ka = [0, tt]. We first observe that dF/d(p = 
at (fi = 0, 7r/2. As -F'(O) = 0, at ip = the function has an 
absolute minimum. The second-order derivative at ip = 
7r/2 is negative, and one can simply prove analytically 
that this point is at least a relative maximum. Numerical 
studies show that it is an absolute maximum, such that 
the smallest transverse frequency is found at wave vector 
fco = tt/o and takes the value 



2 



(24) 



This identifies the soft mode. The corresponding eigen- 
modes exhibit a periodic deformation of the chain at peri- 
odicity 2a, analogous to the zigzag structure. We denote 
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FIG. 7: (color online) Branches of the excitation spectrum 
at (a.) Vt = + 0"*" (just above the critical value) and (b.) 
i/t ~ i^t + 0~ (just below the critical value). In (b.) the 
equilibrium structure is a zigzag, the periodicity is doubled 
with respect to the linear chain and the new Brillouin zone 
is halfed. The four branches of the spectrum are obtained at 
this point by "folding" the two branches of the linear chain 
in (a.). The units and style codings in (a.) and (b.) are the 
same as in Figs. |4] and [G] respectively. Here, i^j'^' ~ 2.05ii;o- 

by bo the amplitude of its oscillations, with bo = 6o(^'t), 
such that the transverse oscillations along y of the ion j 
are described by the function 

ysoft ^ (_i)^6„/2. (25) 

In the following we assume zero temperature and study 
the equilibrium position of the crystal with the transverse 
frequency varying in the interval [vf'^'^ — bv^vf^ + &v\^ 
thus on both sides of the critical point, whereby bv is 
a small positive quantity. Following Landau theory, we 
demonstrate that the zigzag mode of the linear chain, 
given by Eq. (P5|) . is indeed the soft mode, driving the 
instability across the critical point, and we evaluate the 
critical exponents for some quantities of interest. 

Equilibrium positions around the critical point 

In order to determine the behaviour at the critical 
point, we first expand Eq. ([2|) till the fourth order around 



the equilibrium positions of the chain, V = X]/=i 
where I labels the order. The zero order term at leading 
order in is (27l | 

^'"' = ^(A^-l)(7-ln2 + ln(7V) + 0(^^)) ' 

(26) 

where 7 = 0.577216 ... is Euler's constant. The first 
order term vanishes as a result of the requirement that 
we are looking for a minimum. Using the decomposition 
into the eigenmodes of the linear chain, Eqs. (fT3|) -(fT5 |l . 
the quadratic term of the expansion of potential (21) (with 
the use of P^)) and piH12|l ) takes the form 

fc>0.s=± 

(27) 

where is given by Eq. pT|) . while 

and it coincides with Lo\(k), Eq. (H^), for vt > vf''. The 
third and fourth order terms, obtained by using this de- 
composition, arc presented in App. . 

The linear chain becomes mechanically unstable when, 
by varying 1/4, the frequency of the mode with wave vec- 
tor fco = ""/fli Eq. (|25p . vanishes. Starting from this ob- 
servation, we study the behaviour of the corresponding 

(c) 

mode close to the instability point, when vt — Vt ■ For 
convenience, we denote by and ^'q the zigzag modes 
of the linear chain along the y and z direction, respec- 
tively, at wave vector /cq. Around the instability point 
these modes will be coupled significantly to other quasi- 
degenerate modes by the third and fourth order terms 
^(3) and V^^^. These quasi-degenerate modes are long 
wavelengths axial modes at wave vectors (5fc, such 
that \bk\a <C 1, and short wavelength transverse modes 
^fco-i-i5fe' wave vector A: = fco -I- with \bk'\a <C 1. 

At first order in the small parameter \bk\a <C 1, the 
part V^^ , that contains the summands of the third order 
term V'^^^ giving the coupling of the mode at fco with the 
other quasi degenerate modes, has the form 



"^^^ = E E *o (^^^u' + ei^^*^W) + OiSk^a^) , (29) 

^ ^ i5fe>0 cr=y.z 

where we adopted for convenience the notation ^P^^, :— "^k^k^^sk- We note that Eq. is of first order in Sk. The 
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part Vjf^'* of the fourth order term y , which is relevant to the dynamics of the soft mode at ko , involves only the 
transverse modes that are close in fc to fco, and has the form 



(4) 



Sk 



Sk 



Sk>0 



a=y,z 



Sk>0 



E *o^E 



Sk 



Sk 



Sk2 ' ^<5fci+5fc5 



0((5fc"a^) ,(30) 



where A is calculated from the coefficients of V^^'' at kj = 
^0 (j = 1,2,3,4), see App. , and it takes the form A = 
A/N with 



A 



--C(5)^ 
2 32''^^^ a5 



(31) 



The function T in Eq. ([30]) contains a sum of products 

of three amplitudes 'f^^^'' for Jfc ^ 0, and it is of no im- 
portance for the following considerations. The numerical 
factors appearing in Eq. (|30p . multiplying each term of 
the sum, account for all possible permutations of the am- 
plitudes ^"^^ in each summand (see App. and Eq. ((T7|) 
therein). The coupling between the transverse modes 
at ko and the axial modes does not appear explicitly in 
Eq. ([30]) . as it scales with {Sk a)^ ^ 1, and it is hence 
of higher order with respect to the coupling among the 
transverse modes. Since the third order term, Eq. ()29p . 
scales with 6k a, at zeroth order in the expansion in 
\Sk\a and close to the instability, the effective potential 
describing the dynamics of the mode at fcg is given by 
(see Eqs. ^ and ^) 



m 

y 
E 

5fc>0 



(32) 



Sk 



Sk 



Vn 



(4) 



<j=y,z 



where Psk = f3{kQ — 6k) . We now allow the transverse fre- 



quency vt to take values in the interval [vj:'^^ —6v, ly]:''' +61^], 
such that (3sk may take on small but negative values. We 
first determine the amplitude of the zigzag mode fco and 

(c) 

then show that in the vicinity of the frequency z/j no 
other modes arc stable. For this purpose for (3sk < 
we determine the corrections ^'^i^'' to the equi- 

librium positions of the linear chain using Eq. p2p . as- 
suming that these give rise to a small displacement b 
with respect to the equilibrium interparticle distance a, 
6 <C a. In particular, following our hypothesis that close 
to the transition point the soft mode is unique, and it is 
the zigzag mode, we consider the set of solutions where 



(±) 

Sk 



for 6k > 0, and introduce the Fourier amplitude 



of the displacement in the transverse plane g = \/Nb/2, 
as indicated from Eq. p^ . such that 



From Eq. 



one finds g = for /?o > 0, while for /5o < 



Q - 



-N 



p \ 1/2 
AA 



(33) 



This is indeed a minimum if we ignore terms in Vcs with 
non-zero 6k. It will be shown in what follows that this 
minimum is stable with respect to addition of such terms. 

We now demonstrate that Eq. p3|) is actually the 
transverse displacement, giving the equilibrium trans- 
verse positions of the zigzag structure, by verifying that 
Eq. ([33]), together with = for 6k > 0, yields a sta- 
ble solution. To check stability the matrix of the second 
derivatives of Voff with respect to the various variables 
should be calculated. The second derivative of Vcs, given 
in Eq. p2[) . with respect to g is positive, 



dg^ 



= -2m/?o > . 



{e,*5)=}={e,o} 



(34) 



In order to investigate the coupling of the soft mode with 
the modes with 6k ^ 0, one can calculate the second 
derivatives of Vcs with respect to ^'g . We find 



dnd'^sk 



= 0. 



{e,*afc}={e,o} 



This result shows that the derivatives with respect to 
form a sub-block of the stability matrix that can be 
diagonalized separately. All its eigenvalues are found to 
be positive. The other second derivatives at these points 
read 



9* 



K±)2 

Sk 



Sk 



{e,3'5fc}={e,o} 



{e:i'5fc}={e,o} 



{e,*5fc}={e,o} 



{e,*sfc}={e,o} 



mPsk + ^Ag^ 
mPsk + SAg^ 







16A^m , 



WA-^^) 
16A^^} 
(37) 
(38) 



where we have used that > (isk > Po- This result shows 
that the modes ^'^{.^^ and ^'^^^^ are coupled in pairs. All 
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contributions resulting of differentiation of the function 
T in Eq. pop vanish. The stability matrix splits into 
2x2 blocks that can be diagonalized separately. Using 
Eq. ([55]) and /Sq < (i&k one finds that the eigenvalues of 
each block are m{Psk — '2Po) and m{(3sk — ^fio), hence both 
positive. Therefore, a gap opens between the soft mode 
frequency and the frequency of the modes at ^ 
in the vicinity of the transition point. Therefore, the 
instability is driven by the soft modes with wave vector 
fco, determining the order of the zigzag phase 34, 35| . 



where C = 112 C(3)/[93 C(5)], and whose second deriva- 
tive with respect to vt is clearly discontinuous at the 
critical point. This result is consistent with the result 
presented in (23l |. where a discontinuity in the second 
derivative of the ground state energy with respect to the 
particles density was found. 



Discussion 



Behaviour at the critical point 

From Eq. ([32|) . using the results of the previous sec- 
tion we can now write the effective potential for the soft 
modes, which reads 

where A is given by Eq. (j3ip and 



,(c)2 



(40) 



Here, we have used that I3q = u}\ ^^^^ which in turn is 
given by Eq. (l24|). Hence, for V > the potential 
has a single minimum with = 0, and the linear chain 
is the ground state structure, while for V < the po- 
tential landscape has the characteristic form of a Mexi- 
can hat with degenerate zigzag ground states at differ- 
ent angles around the symmetry axis. Indeed, while the 
order parameter q is fixed by condition ()33p . the ratio 
^o/^o is arbitrary. The system hence possesses "Gold- 
stone modes" at zero frequency, which are a consequence 
of the symmetry by rotations around the trap axis. 

The transverse displacement from the trap is given 
from Eq. ((331) by using Eq. ([25]) and the relation 6/2 
g/^fN , which links the displacement in real space with 

(c) 

its Fourier decomposition. Hence, for vt < the trans- 
verse displacement from the trap depends on vt as 



(41) 



with h = 2mv[^^ / A. This behaviour is in agreement 
with the numerical results in [2^ . 

From Eq. (|4T|) we evaluate the difference between the 
ground state energy of the linear and of the zigzag struc- 
ture. Considering the energy per particle, from Eq. (j39p 
we find 



AE 



N 



-mCa {vt - vl ') 



(42) 



Using symmetry arguments we have demonstrated that 
the transition from a linear chain to a zigzag structure, in 
a system of anisotropically confined charges, is a second- 
order phase transition, whose order parameter is the dis- 
placement from the trap axis. This theory has been de- 
veloped in the thermodynamic limit, fixing the interpar- 
ticle distance a as the number of ions was let to infin- 
ity. In this limit, we found that the soft modes are the 
zigzag modes of the linear chain, whose periodicity is 
equal to twice the interparticle distance a. The insta- 
bility is thus driven by these modes as the transverse 
potential is changed across the critical value . 

These considerations are strictly valid for iV cx), but 
can still be useful for finite systems, and in particular 
when the ions are confined in a trap, which provides also 
axial harmonic confinement. While detailed quantitative 
predictions can be only made by accurately evaluating 
the finite-size corrections, we can still make some reason- 
able conjectures, based on previous results in the litera- 
ture and on our theory. Inside a harmonic trapping po- 
tential, the interparticle distance between the ions varies 
along the chain and it is minimal at the center. Numer- 
ical results, based on molecular dynamics simulations, 
showed that in this case the zigzag structure appears at 
the center of the chain where the density is highest 22 1 . 
Analytical studies found that the short wavelength modes 
are characterized by largest displacements at the chain 
center, while the ions at the edge almost do not move [2^. 
In this case, hence, we can still identify the zigzag mode 
of the ion chain with the soft mode. In the presence of 
axial confinement, however, both the transverse as well 
as the axial equilibrium points will change. In particular, 
when going to the zigzag structure the axial density of 
ions in the center will increase. Close to the transition 
point, one finds that the axial corrections to the linear 
chain positions are much smaller than the transverse dis- 
placements from the trap axis as this is a quantity that 
follows the order parameter. This is also confirmed by the 
analysis made for the simple case of three ions in Sec. , 
where close to the critical value of the aspect ratio the 
transverse displacement varies faster than the axial one. 
Therefore, we expect that our theory will still provide 
reasonable predictions close to the critical point, also in 
presence of axial confinement. 
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v^t"'^ (MHz) 




10 10 

a {fim) 

FIG. 8: Phase diagram close to the hnear-zigzag transition 
in the thermodynamic limit, for '"'Ca"'" ions. The horizontal 
axis is the interparticle spacing a in ^m, the vertical axis the 
corresponding critical frequency vj:'^^ in MHz. This graphic 
does not report further curves in the left region, giving the 
transition to more complex structures. A detailed study of 
the transitions to these structures can be found in p^ . [20| ] 
and fH. 
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EXPANSION ABOUT THE EQUILIBRIUM 
POSISIONS OF THE LINEAR CHAIN 



In this Appendix we evaluate the higher order terms 
of the expansion of the potential in Eq. ^ about the 
equilibrium position of the linear chain. For this pur- 
pose we rewrite the interparticle distance as jrj — rj| = 



Hj r Tij + with 



CONCLUSIONS AND OUTLOOK 

The structural phase transition from a linear chain to 
a zigzag configuration, in a system composed of trapped 
singly-charged particles, is a second order phase tran- 
sition. Using a mean field approach we have derived 
a classical model, describing the system at the critical 
point and its vicinity. Our theory is analytical and its 
predictions agree with the numerical simulations of 
and 0. 

The corresponding phase diagram is shown in Fig. [H 
it shows the regions of stability of the linear chain as a 
function of the interparticle spacing a and the transverse 
frequency Vf The phase diagram is evaluated in the ther- 
modynamic limit, corresponding in keeping a fixed as N 
and the chain length go to infinity. The analysis is valid 
for T = 0, where long-range order in one-dimensional 
structures exists. The quantum statistics of the particles 
at these densities seem irrelevant even at these ultralow 
temperatures since the interaction energy at all stages 
much larger than the kinetic energy, and the particles can 
be considered distinguishable at all effects [36|. On the 
other hand, at the critical point, where large fluctuations 
of the transverse motion classically occur, quantum ef- 
fects may be relevant and could be in principle observed. 
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= (<?,: - +Sfj , 

and Sfj — [ui ~ yj)'^ + {zi — ZjY and we have used Xj = 
ja + qj. We now expand in the parameters ey and Tij, 
assuming that they are small with respect to Aij, i.e., 
to the axial equilibrium distances between the ions when 
the chain is stable. We will check later for consistency of 
this assumption. We can write 



(43) 



with 



(0) 



(1) _ 



(i - j)2a2 



(2) _ 



(2(g,:-g,)2-4), 



(3) 



— (<Z.-g,)[34-2(9.-g,f 



\^{lst, + {q.-q,r-3Sl{q,-q,f), 



\i-j\^a^\8 



where we have introduced aij = (i — — j\- We notice 
that 



1/(1) = - 



= 



(44) 



as one can easily verify by using the definition of aij. 
This is satisfied also in the ion chain in presence of an 
axial trapping potential, since V^^^ = determines the 
equilibrium positions. 
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POTENTIAL TERMS FOR THE NORMAL 
MODES OF THE LINEAR CHAIN 

Here we report the third and fourth order terms of the 
potential, decomposed into the modes &k = ©i;"''''"*©!;"'': 
= ^'^(+' - i^'^^"''. The third order term takes the 
form 



= J2 B{h,k2,k3) 

fel+fe2 + fe3=0 



(45) 



where the sum runs over positive and negative values of 
kj and 



B{ki,k2,k3) 



m>0 p— 1 



Term (|45p is real, as it is visible by using the decom- 
position into even and odd modes. In particular, it has 
odd parity, coupling either three odd modes or two odd 
modes with an even one. 
The quartic term reads 



(47) 



^ ^ A(fci,fc2,fc3,fc4) 

fcl+fc2 + fe3+*;4=0 



3 E *fe,*fc.©fc3 0fc4 



with 



A(fci,fc2,fc3,M=-^Ein 



£Q2 

N a' 



sm 



kpTTia , 



•nv' 

7n>0 P— 1 



-(48) 



This term is even, and it thus couples either four modes 
with the same parity, or two odd modes with two even 
ones. 



[1] D. H. E. Dubin and T. M. O'Neil, Rev. Mod. Phys. 71, 
87 (1999). 

[2] F. Diedrich, E. Peik, J. M. Chen, W. Quint, and 

H. Walther, Phys. Rev. Lett. 59, 2931 (1987). 
[3] R. Blumel, J. M. Chen, E. Peik, W. Quint, W. Schleich, 

Y. R. Shen, and H. Wahher, Nature 334, 309 (1988). 
[4] R. Blumel, C. Kappler, W. Quint, and H. Walther, Phys. 

Rev. A 40, 808 (1989). 
[5] J. N. Tan, J. J. Bollinger, B. Jelenkovic, and D. J. 

Wineland, Phys. Rev. Lett. 75, 4198 (1995). 



[6] W. M. Itano, J. J. Bollinger, J. N. Tan, B. Jelenkovic, 
X. P. Huang, and D. J. Wineland, Science 279, 686 
(1998). 

[7] L. Hornekffir, N. Kjargaard, A. M. Thommesen, and 
M. Drewsen, Phys. Rev. Lett. 86, 1994 (2001). 

[8] N. Kjfflrgaard and M. Drewsen, Phys. Rev. Lett. 91, 
095002 (2003). 

[9] A. Mortensen, E. Nielsen, T. Matthey, and M. Drewsen, 
Phys. Rev. Lett. 96, 103001 (pages 4) (2006). 

[10] J. 1. Cirac and P. ZoUer, Phys. Rev. Lett. 74, 4091 (1995). 

[11] F. Schmidt-Kaler, H. Haffner, M. Riebe, S. Guide, 
G. P. T. Lancaster, T. Deuschle, C. Becher, C. F. Rods, 
J. Eschner, and R. Blatt, Nature 422, 408 (2003). 

[12] D. Leibfried, B. DeMarco, V. Meyer, D. Lucas, M. Bar- 
rett, J. Britten, W. M. Itano, B. Jelenkovic, C. Langer, 
T. Rosenband, et al.. Nature 422, 412 (2003). 

[13] C. Wunderlich, in Laser Physics at the Limit (Springer, 
Heidelberg, Berlin, New York, 2001), p. 26, URL 
eirXiv . org : quant-ph/011 1158 

[14] D. Porras and X FCirac, Phys. Rev. Lett. 92, 207901 
(2004). 

[15] M. Pons, V. Ahufinger, C. Wunderlich, A. Sanpera, 

S. Braungardt, A. Sen(De), U. Sen, and M. Lewenstein, 

Phys. Rev. Lett. 98, 023003 (2007). 
[16] J. Taylor and T. Calarco, arXiv:0706.1951 (2007). 
[17] T. Dauxois, S. Ruffe, E. Arimondo, and M. Wilkens, 

Dynamics and Thermodynamics of Systems with Long- 

Range Interactions, Lecture Notes in Physics 602 

(Springer Verlag, New York, 2002). 
[18] R. H. Baughman, S. O. Dantas, S. Stafstrom, A. A. Za- 

khidov, T. B. Mitchell, and D. H. E. Dubin, Science 288, 

2018 (2000). 

[19] G. Birkl, S. Kassner, and H. Wahher, Nature 357, 310 
(1992). 

[20] I. Waki, S. Kassner, G. Birkl, and H. Walther, Phys. Rev. 

Lett. 68, 2007 (1992). 
[21] T. B. Mitchell, J. J. Bollinger, W. M. Itano, and D. H. E. 

Dubin, Phys. Rev. Lett. 87, 183001 (2001). 
[22] J. P. Schiffer, Phys. Rev. Lett. 70, 818 (1993). 
[23] G. Piacente, I. V. Schweigert, J. J. Betouras, and F. M. 

Peelers, Phys. Rev. B 69, 045324 (pages 17) (2004). 
[24] L. D. Landau and E. M. Lifshitz, Statistical Physics 

(Pergamon, Oxford, 1958). 
[25] A. Steane, Appl. Phys. B 64, 623 (1997). 
[26] D. F. V. James, Appl. Phys. B 66, 181 (1998). 
[27] D. H. E. Dubin, Phys. Rev. E 55, 4017 (1997). 
[28] G. Morigi and S. Fishman, Phys. Rev. Lett. 93, 170602 

(2004). 

[29] G. Morigi and S. Fishman, Phys. Rev. E 70, 066141 
(2004). 

[30] D. H. E. Dubin, Phys. Rev. Lett. 71, 2753 (1993). 

[31] G. Morigi and S. Fishman, Journal of Physics B: Atomic, 

Molecular and Optical Physics 39, S221 (2006). 
[32] M. G. Raizen, J. M. Gilligan, J. C. Bergquist, W. M. 

Itano, and D. J. Wineland, Phys. Rev. A 45, 6493 (1992). 
[33] N. Ashcroft, Solid State Physics (Holt-Saunders Int. 

Eds., Tokyo, 1976). 
[34] O. Biham, D. Mukamel, and S. Shtrikman, Phys. Rev. 

Lett. 56, 2191 (1986). 
[35] D. Mukamel and M. B. Walker, Phys. Rev. Lett. 58, 2559 

(1987). 

[36] J. Yin and J. Javanainen, Phys. Rev. A 51, 3959 (1995). 



